On a certain analogy between hydrodynamic flow in porous media and heat 

conductance in solids 
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We consider a porous medium being saturated with a pore fluid (Biot's theory). The fluid 
is assumed as incompressible. It is shown that the general integral of the elastic and pressure 
equations can be written in form of a time dependent vectorpotential F being a solution of 
a homogeneous, fourth order differential equation. The obtained equation for F is of a more 
general form than the corresponding thermo-elastic vectorpotential, being a solution of a time 
dependent and inhomogeneous vector bi-Laplacian. Both vectorpotentials do, however, agree 
for stationary problems in general and for certain particular boundary conditions (irrotational 
deformations). An example of an irrotational deformation is studied in detail, exhibiting known 
properties of classical vector diffusion. 

PACS number(s):?? 



o 

CO 



I. INTRODUCTION 
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Imposing an inhomogeneous pressure/temperature 
field on a (porous) material causes in general mechani- 
cal deformations and stresses. These 'induced' elastic 
deformations in turn can and do influence the fur- 
ther mass/heat transport. The problem is described 
as a coupled field problem consisting of (at least) two 
equations (a) an clastic equation describing mechan- 
ical displacements due to a given source field (pres- 
sure/temperature gradient field) and (b) at least one 
transport equation for mass/heat with a source term 
originating from the mechanical deformation (diver- 
gence of the displacement field). 

The entire description bases on continuity equations 
and constitutive laws which are discussed in textbooks 
in great detail IjJ]. Perhaps less known is Biot's theory 
of hydrodynamic flow in porous elastic media . We 
will address this problem below in more detail. In par- 
ticular we wish to draw some attention on similarities 
for the flow and the heat problem. 



II. GENERAL CONSIDERATIONS 

We give the basic outline of Biot's theory for the 
particular case of an incompressible pore fluid. This 
considerably simplifies the governing equations and 
their discussion, because all material derivaties can 
be eliminated. First we discuss the basic elastic equa- 
tions. This is followed by a detailed discussion of the 
flow equation for incompressible pore fluids. There- 
after we discuss some properties of the coupled equa- 
tions, and finally establish their general integral fol- 
lowing the spirit of Galerkin's method. We specialize 
the obtained general integral for the case of irrota- 
tional displacement fields and make the connection to 
the problem of thermoelasticity. 



A. Preliminaries 

In the following bold-face letters always indicate 
vectors/tensors. 



The continuity equation for momentum density, 

div a + f = p s ii, (2-1) 

represents the mechanical equilibrium condition, 
where a — (<7;j) is the stress tensor, f any external 
body force, p s the solid's mass density and u = (itj) 
the mechanical displacement. We will neglect the in- 
ertial terms in the following and will consider Biot's 
theory of porous media for the case of symmetric elas- 
ticity, i.e. cry = aji. Deformations are characterized 
by the strains 2ity = (dui/dxj+duj jdui). Employing 
Hooke's linear relation between stresses and strains 
for isotropic materials one obtains Lame's equation 
for the displacements, 

2 

(1 - 2^)A u + grad div u + — (1 + v)(l - 1v) f = 0, 

E 

(2.2) 

with Poisson contraction number v and Young modu- 
lus E. There are a few quite general results about this 
equation, which will be useful in the following. Taking 
the divergence of Eq.(2.2) one has, 



A div u 



(1 - v)E 



div f , 



(2.3) 



i.e., div u is the solution of a Poisson equation and for 
vanishing body forces f it is even a potential function. 
The physical significance of div u is that it measures 
the local relative volume change due to the deforma- 
tion which is an invariant. Pure shear deformations 
are unreflected in Eq.( |2.3| ) because they preserve the 
volume. Of course all information on vorticity has 
been lost by passing from Eq.(2.2) to Eq. (|2.3|) . On 
the other ha nd t aking the Laplacian on Eq.( |2.2| ) and 
inserting Eq^^) into yields, 



A A u = 



2(1 + 1/) 
E 



1 



2(1-1/) 



grad div f — A f 



(2.4) 
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This equation is of fourth order while Eq.(2.2) is of 
second order; hence their solutions will be in general 
not identical. In his seminal work Galerkin resolved 
this problem. He introduced another vector function 
F satisfying the (inhomogeneous) biharmonic equa- 
tion, 



(2.5) 



A AF = --{l + y) 



with f being the same body force as in Eq. fl2.2p . He 
showed that the integral of Eq.(2.2) can be expressed 
in form of derivatives of F, 



A F 



2(1-1/) 



grad div F. 



(2.6) 



We would like to note that the foregoing consider- 
ations made no particular assumptions on the body 
force field f in particular they apply to conservative 
fields. Typically f is explicitly known like for gravita- 
tional forces, or it is itself a solution of another inde- 
pendent equation. The latter is the case for thermo- 
elastic problems, in which 



f = —a K grad T 



(2.7) 



and 



J^(</W ) + ( <7'/v j = i). 



(2.11) 



respectively, where p s and pi denote the solid and fluid 
mass densities. In the following we make the two im- 
portant assumptions that the mass densities are con- 
stants (p s = cons t, and pi = const.). 



From Eqs.fl2.9j), fl2.10| ) and fl2.1l|) and the assump- 
tion of constant mass densities follows a) that the local 
volume is conserved because the inflow of fluid volume 
into a volume element equals the solid volume outflow, 
and b) that the flow equation takes the following sim- 
ple form (n/p — const), 



— A p — div u, 



(2.12) 



i.e. a Poisson equation for the pressure with the rate 
of the relative elastic volume change as source term. 
The porosity <fi does not enter explicitly into the flow 
equation because of constant mass densities. As soon 
as u has b een determined from the coup led E qs.(2.2) 
and ( 2.12 ) </> can be computed from Eq.( 2.10 ). 



where a is the heat expansion coefficient, the K = 
3 ^ 1 E 2l/ ) accounts for the material's compressibility, 
and T is the temperature field. Typically Eqs. (|2.2|) 
and (2.7) are supplemented by a heat conductance 



equation, 



(2.8) 



with x being the coefficient of heat conductance di- 
vided by the specific heat capacity. Equation (2.5) 
then becomes a biharmonic equation with time de- 
pendent -but known inhomogeneitis- which do not 
dependent on the deformation themselves. As will be 
shown further below does flow in porous materials im- 
ply in general 'body forces' being deformation depen- 
dent themselves. The deformation acts back on the 
flow. We will come back to this point in more detail 
further below. 



2. Flow equation 

Lets consider fluid flow in a porous continuum. The 
porosity (local volumetric fluid fraction) can be de- 
scribed by a scalar field </>(r, t) . It is common to as- 
sume that the relative fluid-solid velocity 0(v — u) can 
be expressed by Darcy's relation, 



f = gradp = <j)— (v — u), 



(2.9) 



where v and u are the fluid and solid velocities, k the 
mechanical permeability, p the fluid viscosity, and p 
the hydrostatic pressure field. 

The continuity equations for solid and fluid mass 
are, 



d (, 

mil 1 



div 



((1 - 4>)Ps<i 



0, (2.10) 



B. The coupled system 

1. The Galerkin- Biot potential 

After the a bov e prep aratio ns we will focus on the 
coupled Eqs.(2.2) and ( 2.12 ) describing elastic-flow 
interactions for incompressible pore fluids. Note that 



grad p, 



(2.13) 



in Eq.(|2_2| 

equation 




comp. Eq.(2.9). In general Lame's vector 
2|) is much harder to solve than the scalar 
However, 



Poisson Equation ( 2.12 ). However, for the very par- 
ticular case of irrotational displacement s Lam e's equa- 
tion can be integrated easily, see Sec. Ill A. In gen- 



eral (and also for most boundary conditions of prac- 
tical interest) displacements and velocities will have a 
non-zero rotational part. Therefore div u or its time 
derivative will not contain the full information about 
displacements and strains, i.e. all elastic volume pre- 
serving deformations (shears) are unreflected in it. 

It has been convincingly demonstrated that for 
heat-conductance in solids the relative volume 
changes, div 
However, 



u, actually do also appear in Eq.(2.8). 
the corresponding coefficient is for most 
solids so small that this term can be safely disregarded 
simplifying the theoretical treatment considerable Q . 
As an analogy to heat-conductance in solids one can 
think of a strongly compressible fluid, i.e. a gas in a 
porous medium. One can neglect the back-influence 
of the pressure induced elastic distortion on the pore 
pressure field and is basically left with the equations 
decribing a corresponding thermo-elastic problem. On 
the other hand does flow of an nearly incompressible 
fluid in a porous medium represent the 'opposite' case: 
in order flow to happen elastic deformation must take 
place. Without elastic deformation there is no flow 
and vice versa. 
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We reconsider the flow problem Eq.(2.5) and make 
some small refinements. Of course the pore pressure 



held is unknown in Eq.(2.5) since we h ave not incorpo- 
rated any information from Eg. ( 2. 12 ). This needs to 
be done in order to obtain a fully dec oupl ed equation 
for displacements. Starting from Eq.( 2.12 ) we have, 



f = — grad p = — grad 

K 



div u(r') 



, 47r|r — r 



jrdV. (2.14) 



with D being already defined in Eq.( 2.16| ) The dis- 
persion relation has the form of a diffusion only prob- 
lem, i.e., all solutions are damped in time and there- 
fore the equations are stable. The damping is con- 
trolled by a constant of diffusion D in a rough esti- 
— i 2 

mate of order 10 — which is orders of magnitude 
larger than for molecular diffusion (we estimated D 



from k 10 15 
v = 0.25). 



E « 10 11 Pa, n w 10~ 3 Pa s and 



This equation still needs to be complemented by 
boundary conditions for the pressure field. Insert- 
ing Eq.(2.6) into Eq.( 2.14 ) and integrating employ- 
ing Helmholtz's theorem Jl] one obtains a relation be- 
tween f and F which is in turn inserted in Eq.(2.5), 



A A F = — grad div F, 
D 



with the abbreviation 



D= K - 



E(l 



(2.15) 



(2.16) 



Equation (2.15) is the main general result in this 
section. It needs to be complemented by boundary 
and initial conditions for displace men ts an d pre ssure. 
It is not difficult to see from Eqs.(2.6) and (2.15) that 
if curl u = the displacements are characterized by a 
vector-diffusion equation. However, the general case 
with non- vanishing curl is much more complicated (it 
is described by an intcgro-differcntial equation in u). 



2. Dispersion relation 

While theoretical simplifications within models are 
usually appreciated a remaining question is whether 
the simplified model assumptions do lead to a consis- 
tent physical picture concerning certain criteria. One 
such criterion is stability. For example, is it obvious 
that the model equations are stable if one disregards 
inertial forces for the elastic continuum phase? This 
problem is most convenient considered investigating 
the dispersion relation w(k) obtained from an expan- 
sion after eigenfunctions. For brevity we only present 
the basic result in the following. 

Assuming certain mathematical integrability condi- 
tions on displacements u(r, t) and pressure p(r, t) both 
fields can be represented as Fourier-Integrals, 



and 



u(r, t) = J dk u(k) e 
p(r,i) = J rfkp(k) e 



i(ujt — kx) 



i(ujt— kx) 



(2.17) 



(2.18) 



Inser ting this two representations into the basic 
Eqs.( |2.2| ) and (2.12) one obtains the dispersion rela- 
tion, 



III. THE POINT SOURCE 
HYDRAULIC LOADING) 

In the following we consider a particular simple case 
in three dimension: at time to the material is loaded 
by a pressure peak S(t)5(r). Dimensionless quantities 
for the length scale, stress and time scale will be useful 
(denoted by a star); 



,1/2 ' 



P 



P_ 

G' 



f* 



,G 



t- 



(3.1) 



with k, 2G = E/(l + v) and fj, being the permeability, 
elastic modulus and fluid viscosity respectively. 



A. Pressures and Displacements 

Because the chosen symmetry of boundary condi- 
tions guarantuees the existence of an irrotational solu- 
tion Lame's equation Eq.(2.2) simplifies considerable, 



div* u* = a* (p* — Pq), 



(3.2) 



where p*$ denotes the equilibrium pressure and 2a* = 
(1 — 2^)/(l — v) a material constant. Note that the 
solutio n fo r u* can be obtained by direct integration 
of Eq.(|J), 

„*2 



u*(r*) = grad* 

47T 



dV*' P } } ~f.° . (3.3) 



Equation (3.2) can be inserted into the pressure equa- 
tion (2.12) giving an (inhomogeneous) heat conduc- 
tance equation, 



A*p* - a* 2 ^ = -S(t*)S(r* 



(3.4) 



with three dimensional solution (Green's function) 
1 



p*(r*, t*)-p* 



* 2 



(A/n) 



3/2 -At" 



where, 



A = 



At* 



(1 - 2v)[i 
8(1 - v)Gt 



(3.5) 



(3.6) 



Thus incompressible flow in a porous medium whose 
displacement field is irrotational corresponds to a heat 
conductance problem in a fiktive non-porous medium. 
One can insert the pressure field Eq. Q3.5D into Eq.(3.3) 
and evaluate the integral, 



1 f 1 f^' 1 "*' 



dz - ^re~ A * 
r* 



Lj(k) = i D k z 



(2.19) 



(3.7) 
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In Fig.|| we show the dimensionless radial displace- 
ments, u*, for two different times. 



0.040 



.. 0.030 - 



S 0.020 - 



Bi 0.010 - 



0.000 




2.0 3.0 
Radius r 



FIG. 1. Radial displacement according to Eq.(3.7) for 
two different times A — 1.0 and A = 0.5. 



FIG. 2. Scaling function H(S) a s fu nc tion of the re- 
duced radius 5 = \f~Ar *, comp. Eqs.Q, and 
The function, being essentialy a rescaled radial deforma- 
tion, shows four regions of interest: (a) for small radii the 
positive large but finite deformation, (b) for So ~ 1 a root, 
(c) a negative, minimum deformation at <5 C as 1.6 and (d) 
a decay towards zero for large radii. 

Let So denote the root of the forgoing equation. We 
find, 
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At 



(3.10) 



within three percent error. The volumes exhibiting 
positive and negative stresses are therefore separated 
by a surface of radius r^ at time A around the center, 
or equivalently tq = 2\J Dt with diffusity D defined in 
Eq. (2.16), eg. an order of magnitude estimate gives 
ro Ri 20 cm for t = 0, 1 sec for the neutral zone. An- 
other interesting property is the asymptotic behaviour 
of the deformation for <5 — * oo. The asymptotics is 
given as, 



The displacements being zero at the symmetry cen- 
ter raise to a maximum value and drop for larg er r adii 
towards zero again. As can be seen from Eq.(3.7) do 
the displacements follow a scaling form. The scaling 
variable i5 = >/Ar* may be anticipated for a diffusion- 
related problem. 



B. Deformations 

It is obvious that the deformations can be written 
in scaling form, 



with H{8) being a scaling function of the form, 



H{6) = e 



1 



~ 6 ' " - / dze- z ' 



5 2 S 3 J 



(3.8) 



(3.9) 



being plotted in Fig. |[ 
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(3.H) 



It is interesting to note that neither elastic constants 
nor fluid properties enter this relation. Only the 
permeability enters this compressional far-field decay. 
This relation should be compared to the simple case 
of a cavity of radius R bearing pressure — p within an 
infinite non-Biot medium with vanishing pressure at 
infinity, 



P_ 

' 2 



i?\3 

r , 



(3.12) 



when one measures the bearing pressure in units of 

r. Both relations are essen- 



the shear modul G = 2 (i+i>) • 
tiall y the same, where the typical length-scale re 1 / 2 in 
Eq.(3.11) given by the solids permeability acts similar 
to the cavity-size in non-Biot elasticity. In principle 
this should allow to measure the permeability exper- 
imentally for S — > oo, which means in practice for 
distances r » 2yDt. This observation leads us to 
the conclusion that the more complicated 'Biot elas- 
ticity' reduces for length scales much larger than the 
scale of diffusion towards convential linear elasticity. 

The maxima of the tensile deformations, S — > 0, are 
readily calculated using 



, du* 
hm — — 

S-vO dr* 



A 3 ^ 2 

3^ 3 / 2 



(3.13) 



This value does not depend on the material's perme- 
ability. 

Before turning towards the mechanical stresses we 
quickly summarize the results for the polar deforma- 
tions u*,,. Fig. U shows the polar deformations, in a 
rescaled form. 
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FIG. 3. Polar deformation = u*/r*. Because the 
materials moves outwards from the injection center all de- 
formations are positive. Note that the polar deformations 
are monotonically decreasing for increasing 8 while the ra- 
dial deformations exhibit a negative minimum.. 



One has, 



1 /A\3/2 / 

2Kn) ( £ 



H(S) 



(3.14) 



with H(S) and 6 as defined in Eq. (3.9). The limiting 
cases are, 



lim liL = lim u* r , 



(3.15) 



and 
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FIG. 4. Rescaled radial stress E* r = 

as a function of the rescaled radius 5 for two values of Pois- 
son numbers. For v — > 1/2 the scaling function decays like 



which is positive only. 



The characteristic points of the rescaled stress E* r 
do now depend on v albeit not very strongly, i.e., 
for lower v one finds lower maximum tensile stresses, 
lower roots So and higher compressive stresses. Gener- 
ally spoken do the polar deformations u*^ strengthen 
the tensile stresses and weaken the compressive region. 
Higher contraction numbers v do stronger weight po- 
lar deformations on radial stresses. The two limiting 
cases for the radial stresses are, 



lim < 

5^0 



2(1 + 1/) /A\3/2 



3(1 



± ,Ay 



(3.19) 



and 



lim u*, 

6^00 ' 



- lim u* 

2 8 — *oo 



(3.16) 



lim a* r 

6 — »oo 



tt(1 - 2v)r* 



(3.20) 



IV. CONCLUSION 



C. Stresses 



We now proceed to calculate the radial stresses a* r 
which we measure in units of the shear modul G, 



( 1 -2v) ~ ^ U * rr + 2 ""^) ' ^■ 17 ^ 



which is equivalent to 
1 -2i//vl\-3/2 



(-) <■=(!- v)H{8) + v{er &2 H{5)). 



(3.18) 



We have plotted the rescaled stresses Eq. (3.18) in 
Fig. [| for Poisson numbers v — 0.2 and v — 0.4. 



We have considered a simplistic model for the fluid 
flow within a porous elastic solid. The presented de- 
scription is given by linear equations. Essential sim- 
plifications are: Vanishing inertial forces of the solid 
phase and vanishing fluid compressibility. The first 
assumption relieves us of the problem of wave propa- 
gation in Biot media. The second assumption allows 
us to construct a linear description. We tried to ob- 
tain some general results for the linear problem. It 
appears that the general integrals for Biot's and the 
thcrmo-elastic problem are different even within a lin- 
ear description. While in both cases the body force 
field is governed by scalar potentials, e.g. pressure 
and temperature respectively, Biot's problem is a fully 
cross-coupled field problem. This changes the struc- 
ture of th e governing Galerkin-Biot vector-potential, 
Eq. (2.15), substantially in such a way that mixed 
space-time derivatives do appear. However, we were 
able to show that an irrotational displacement field 
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garantuees an equivalence of both problems. Though 
irrotational displacement fields are rare in geometries 
of practical interest, we studied the point symmet- 
ric case in three dimensions in order to exemplify the 
fluids action on the elastic strain/stress distribution. 
We quantitatively discussed the situation that arises 
when a pressure pulse within a 'borholc' drives fluid 
through the solids pore-space. We found that radial 
tensile stresses arise in a region where non-permeable 
solids exhibit only compressive stresses. This surely is 
of interest to selected problem in soil meachanics and 
fracture mechanics, i.e. hydraulic fracturing. 
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